Stability and structure of two coupled boson systems in an external field 
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The lowest adiabatic potential expressed in hyperspherical coordinates is estimated for two boson 
systems in an external harmonic trap. Corresponding conditions for stability are investigated and 
the related structures are extracted for zero-range interactions. Strong repulsion between non- 
identical particles leads to two new features, respectively when identical particles attract or repel 
each other. For repulsion new stable structures arise with displaced center of masses. For attraction 
the mean-field stability region is restricted due to motion of the center of masses. 
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Introduction. Bose-Einstein condensation is by now 
routinely realized in many laboratories P, Q- Static 
and dynamic properties are currently investigated by use 
of Feshbach resonances to manipulate the effective two- 
body interactions. Two-component and mixed systems 
open for a variety of combinations. Two-component con- 
densed bosons systems consisting of the same atoms in 
different spin states are formed (41 and collective oscil- 
lations studied for ®'^Rb @. Condensates of two different 
species were recently created combining ^^Rb and ''^K ||6j. 
Also mixed systems of condensed boson and degenerate 
fermions was recently experimentally obtained 7] . 

Theoretical descriptions often use the mean-field ap- 
proximation with zero range interactions J^, _2] . Both sta- 
bility and average properties are successfully described 
for one-component boson systems which so far are the 
most intensively studied. The simplest model for two- 
component and mixed systems is a coupled mean-field 
approximation 0, ^ . Two-cluster features may then 
be accounted for by using mean-field intrinsic and relative 
coordinates. However, for such systems a more compli- 
cated cluster structure may be favorable corresponding 
to lower energy. Competing or maybe co-existing new 
types of structures may appear. The stability conditions 
may be altered when degrees of freedom beyond the mean 
field are allowed. 

The mean-field models are unable to investigate ef- 
fects of correlations Correlations can be studied in 
hyperspherical models, which have been applied to one- 
component symmetric boson systems in external fields 
|12|. An extension to two-component systems requires 
at least introduction of a relative and a size coordi- 
nate for each subsystem. Compared to a one-component 
system the difficulties increase substantially. Therefore, 
before embarking on such a large project, it is advis- 
able and illuminating to study the decisive large-distance 
limit with zero-range interactions and simple approxi- 
mate wave functions. The purpose of the present letter is 
to provide estimates for conditions of stability and occur- 
rence of new structures for two coupled boson systems. 

Theoretical method. The N bosons have masses rui 
and coordinates f^. We assume a division into two sys- 
tems, A and B, of identical bosons with particle num- 
bers Na and Nb {N — Na + Nb) and total masses 



Ma = Natua and Mb = Nbitlb {M = Ma + Mb). 
We shall use hyperspherical coordinates with the hyper- 
radius defined by jlJb il^ 
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where m is an arbitrary normalization mass, = ri—rj 
and R = 'Y^^rriiri/M is the center-of-mass coordinate. 
The remaining degrees of freedom are described by the 
hyperangles collectively denoted by Jl. A system divided 
into two different subsystems of non- identical particles re- 
quires independent parameters for each component. The 
size coordinates defined as the corresponding two hyper- 
radii are then the natural choice analogous to mean-field 
coordinates for each component. In addition we also al- 
low the correlations corresponding to relative motion of 
the two center of masses. This is the minimum number of 
parameters for a two-component system where the sizes 
and center of mass positions can be independently var- 
ied. Thus we define the hyperradii and center-of-mass 
coordinates for the two systems A and B 
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The average separation of the two components is then 
given by the coordinate r = Ra — Rb- The three length 
coordinates, pA, Pb, and r, are related to p by 



mp'^ — mAp\ + mBP%+mrr'^ , urir 
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We aim at using the hyperspherical adiabatic expan- 
sion method. We only include the first adiabatic angular 
wave function, i.e. 
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where the separable center-of-mass motion is removed. 
We first assume that s-waves dominate in the angular 
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expansion and second that $ is independent of the an- 
gles. Then the wave function is automaticaUy symmetric 
under permutations of particles within each subsystem. 
These assumptions are severe, but not unreasonable for 
spatially extended dilute systems. 

The three short-range interactions are all chosen as 6- 
functions with strengths expressed in terms of the corre- 
sponding two-body scattering lengths a a, as, and qab, 
i.e. 
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where fj,AB = rnAniB/i'tnA + ttib)- Integrating these 
interactions over all hyperangles we arrive at the 
Schrodinger equation which determines / in Eq. lO, i.e. 
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where the adiabatic effective potential for Na ^ 1 and 
iVs > 1 is 
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where we have defined the trap length for the external- 
field frequency uj as = h/iniA'-o)- The dimensionless 
function / depends on Na^ Nb, Pa/t, and pb/t. It arises 
from the angular average of the interaction Vab between 
the particles in systems A and B. 

Effective potential. The effective potential Ueff in 
Eq. lfTn|) consists of three types of terms, i.e. from 
the external field (first three), the generalized centrifu- 
gal barrier (next two), and the angular average of the 
interactions and kinetic energies (last three). In total, 
Ueff/(huj) depends on the particle numbers Na and Nb, 
the ratio mA/mB, three coordinates and three scattering 
lengths measured in units of bf. 

The function / in Eq. (|10() is a double integral with 
rather simple properties, i.e. it is monotonous as function 
of r, it vanishes when r > pA + Pb, has a positive max- 
imum and vanishing first derivative at r = 0, and it is 
unity when Na = Nb and pA = Pb- 

The behavior of Ueff is qualitatively different for at- 
tractive and repulsive two-body interactions. The inter- 
action terms {p~^) dominate for small distances, and for 



strong attraction Ueff has no minimum. The system is 
always spatially confined by the external harmonic field. 
When subsystems A and B are identical, the effective po- 
tential has a local minimum when N\a\/bt < 0.67 0,11. 
For asymmetric systems the corresponding stability con- 
ditions are more complicated, as we shall see. 

To illustrate, we show in Fig.^for three repulsive scat- 
tering lengths the minimum of Ueff with respect to pA 
and as a function of r and various divisions of the total 
particle number 40000. For symmetric divisions coincid- 
ing centers of mass (r — 0) are preferred. As the asym- 
metry increases {Nb < OASANa), a deeper minimum 
develops at finite r-values comparable to the trap length. 
The minimum at r = remains stable because the barrier 
separating the minima is much higher than the zero-point 
motion in both wells. The deepest minimum occurs for 
the asymmetric system with Nb/Na ~ 1/40. 
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FIG. 1: The effective potential in Eq. IIOII as a function of 
r/bt for *'^Rb masses niA = ms = 1.44 x lO'^^kg, trap length 
bt = 2.22/^m {lo = 2n X 23.5 Hz), particle numbers A'' — 
Na + Nb — 40000, and oa = as = uab = 103ao, where ao 
is the Bohr radius. The curves correspond to the different 
values of {Na, Nb) given on the figure. 

Structure. The basic characteristics of a system are 
size and energy. The root-mean-square distance is 
the measure of the size of system A defined by 
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where () denotes an expectation value. Correspondingly 
we define mean distances r and fs for the total and the 
B-system, i.e. iVf^ = (p^) and iVsr^ = {p%). 

When the different particles attract each other, i.e. 
O'AB < 0, the effective potential only has the minimum 
for coinciding centers where r = 0. The structures are 
then similar to those of one-component systems. There- 
fore, we shall here concentrate on the different structures 
appearing for repulsive cases where uab > 0. 

The average sizes corresponding to Fig. ^ are shown in 
Fig. For the minimum at r = we obtain large sizes, 
i.e. « fs w f « 3bt. This structure again closely 
resembles the one-component system. For a given divi- 
sion Nb/Na both subsystems contract and stabilize on 
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a lower level as r increases to a few time times the trap 
length. The smaller the system, the larger the contrac- 
tion. 

The effect of the repulsion is that the two systems try 
to avoid each other while staying at relatively small dis- 
tances as dictated by the external field. The large subsys- 
tem {A) is exploiting the space almost like it was alone in 
the trap. When one subsystem [B) is small, it becomes 
advantageous to place about half of its particles outside 
the other subsystem. Then the repulsion on B from the 
trap and the subsystem A is minimized. 
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FIG. 2: The root-mean-square radii ta and rs as a function 
of r for the systems specified in Fig. 

The minima of the effective potentials define the pre- 
ferred structure. As Nb/Na changes, the system with 
the largest particle number {Na) always has the same 
size of about 3bt in these minima. In contrast, the size 
of the other subsystem decreases after Nb < 0.134A'^ 
where the lower minimum at finite r appears. The val- 
ues of fs converge towards bt when Nb/Na 0. 

The center of mass of the small system (B) remains 
within the radius ta until iV^ k, 0.004A'^. For even 
smaller subsystems the center of B moves a little outside 
fyi, but still the distance between the centers at the min- 
imum remains smaller than f^+^s- Therefore, complete 
separation between the systems does not occur. 

Stability. The stability of the system is guaranteed 
for repulsion within both subsystems, i.e. when qa > 
and ub > 0. Instability arises for attraction where one 
or both subsystems can collapse into high-density states 
of lower energy. In Fig. |31 we illustrate these stability 
regions for a specific system as a function of the scat- 
tering lengths a A =0-3 < and uab- The results are 
extracted from the effective potential in Eq. IjlOl) . Sta- 
bility is defined as the existence of a local minimum. We 
first consider r — and search for minima by variation 
of PA and pb- No stationary point exists for A^^ = Nb 
when 
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which follows the critical stability line in Fig. for neg- 
ative flyiB-values. 

Maintaining r = we find that any stationary point of 
Uef f has negative curvature (second derivative) in one of 
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where Eq. lO for < a^s/fot < 0.0001 follows the criti- 
cal stability curve in Fig.|3| Eqs. (221j CSl, and 121) are 

comparable with the mean field approximation 
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FIG. 3: Stability regions for a two-component system with 
Nb = Na = 2000 as a function of ha = as < and aAB- The 
thick solid line separates stable (right side) from unstable (left 
side) regions when we allow three degrees of freedom, i.e. pA, 
Pb and r. The different curves are related to Eqs. 11211 (dot- 
dashed), llOll and |[nt (long-dashed), llT^ and I^SJ (short- 
dashed), respectively. 

When 0.0001 < aAB/bt, the minimum of Uefj occurs 
for finite r, i.e. separation of the center of masses is 
favorable. Furthermore, the curvature at r = in the 
r-direction is negative when 
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where Eq. (|15() is decisive by defining the critical stability 
curve in Fig. E for 0.0001 < aAB/bt < 0.00028. Eq. ^ 
resembles the mean-field expression derived in Q . 

The critical stability curve is in agreement with mean- 
field computations for aAB/bt < 0.0001, but strongly de- 
viating for 0.0001 < aAs/bt where the system exploits 
the r-degree of freedom and becomes unstable for a range 
of attractive a^/^t-values. This is illustrated in Fig.^for 
a set of scattering lengths where a very shallow local min- 
imum is present for coinciding centers, i.e. r = 0. Along 
the line defined by p^i = ps a maximum is found around 
20bt and a minimum at 48.46t. This degree of freedom is 
analogous to the one-component mean-field coordinate. 
The two-component mean-field coordinates are presum- 
ably similar to the use of pA and pb for r = 0. A more 
detailed comparison is required to distinguish. 
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Moving now in the direction where the sum of pA and 
Pb is constant and the difference varies, we find rather 
smah barriers for [pa.Pb) ~ (406t, 50&t), (50&f , 406t). 
This degree of freedom is analogous to the two- 
component mean-field coordinates. We now decrease the 
repulsive scattering length aAB by a factor of two from 
the value in Fig.^while maintaining ua- Then the shal- 
low minimum disappears in agreement with the Fig. |31 
where the decrease of uab corresponds to crossing the 
curve described in Eq. where one curvature changes 
sign. The larger repulsion in Fig. ^stabilizes the system 
as predicted in the mean-field approximation. However, 
if the r-degree of freedom also is allowed, the system col- 
lapses by moving the centers apart from each other as 
described by Eq. (|15|l . The local minimum in Fig. 0] is 
in fact unstable. Stability is regained for a finite r-value 
when aA is less attractive and uab > 0.000286i, see the 
upper middle part of Fig. 
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FIG. 4: Contour plot of the potential Ueff/{Nhij) in Eq. l(T?1jl 
as a function of pA and ps for r = 0, A''^ ~ Nb = N/2 — 
2000, aAB = 0.00026t, and ua ^ ob = -0.00042bt. 

The critical line for stability of each of the inde- 
pendent systems A and B is also shown in Fig. El as 
aA/bt = -0.67/Na = -0.00034. The repulsion between 
A and B stabilizes the total system even when each com- 
ponent would collapse when left alone. The reason is 
that the contraction of each subsystem increases the mu- 
tual repulsion. The balance of these forces results in 



increased stability of the total system for moderate re- 
pulsion. In contrast, less stability results for strong re- 
pulsion where the collapse proceeds by circumventing the 
individual barriers via relative center-of-mass motion. 

The two minima in Fig. ^ are not found when a a = 
aB and Na = Nb- This feature only seems to occur for 
asymmetric systems. When both minima are present, 
each corresponds to a stable configuration. The vibra- 
tional zero-point energy is much smaller than the hardly 
visible, but stabilizing barrier in Fig. ^ Furthermore, 
estimates of the WKB tunneling probability show decay 
rates much smaller than the frequency of the external 
harmonic field. Co-existing structures are then possible. 

Conclusion. We have investigated coupled two- 
component boson systems in a confining trap. The de- 
grees of freedom are the sizes of the two components and 
the centcr-of-mass distance. Two new features emerge. 
First, new stable structures appear in the trap when iden- 
tical and different pairs of bosons all repel each other and 
one of the components at the same time has many more 
particles than the other. Then a finite distance between 
the centers of the two components is energetically favored 
and two stable co-existing structures are possible. 

The second new feature appears when identical bosons 
attract each other and different bosons repel each other. 
Some structures previously computed to be stable in the 
mean-field approximation are unstable. In the present 
model with more degrees of freedom, when the repulsion 
is sufficiently strong, the mean-field minimum becomes 
unstable against relative motion of the distance between 
the two centers of masses. The mean-field regions of sta- 
bility are then further restricted. 

The inaccuracies in the computations are all related 
to the choice of parameter space. Improving by includ- 
ing two-body correlations is possible for example by in- 
creasing the variational space of the wave function as 
done in |llL for a one-component system. These im- 
provements are then correspondingly expected to be sub- 
stantial for hyperradii down to the range of the interac- 
tion, relatively small for hyperradii corresponding to con- 
densate distances, and insignificant for larger distances. 
Thus the qualitative structure would remain unchanged 
and at the crucial distances for condensate formation the 
results are rather close to being quantitatively correct. 
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